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ABSTRACT 



We present a study of the three forbidden oxygen lines [OI] located in the optical region (i.e., 5577.339 A, 6300.304 A and 6363.776 A) 
in order to better understand the production of these atoms in cometary atmospheres. The analysis is based on 48 high-resolution and 
high signal-to-noise spectra collected with UVES at the ESO VLT between 2003 and 201 1 referring to 12 comets of different origins 
observed at various heliocentric distances. The flux ratios Iaoolhm an d GIK=I lsll l{l bi0Q + few) are evaluated to determine the parent 
species of the oxygen atoms by comparison with theoretical models. This analysis confirms that, at around 1 AU, H 2 is the main 
parent molecule producing oxygen atoms. At heliocentric distances > 2.5 AU, the G/R ratio is changing rapidly, an indication that 
other molecules are starting to contribute. CO and C0 2 , the most abundant species after H 2 in the coma, are good candidates and 
the G/R ratio is used to estimate their abundances. We found that the CO2 abundance relative to ILO in comet Q4 (NEAT) observed 
at 4 AU could be as high as ~70%. The intrinsic widths of the oxygen lines were also measured. The green line is in average about 
1 km s -1 broader than the red lines while the theory predicts the red lines should be broader. But at 4 AU, we found that the width of 
the green and red lines in comet Q4 (NEAT) are the same which could be explained if C0 2 is the main contributor. 

Key words. Comets: general - Techniques: spectroscopic - Line: formation 
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1. Introduction 

Comets are small bodies formed at the birth of the Solar sys- 
tem 4.6 billion years ago. Since they did not evolve much, they 
are potential witnesses of the physical and chemical processes 
at play at the beginning of our Solar system (Ehrenfreund & 
Charnley 2000). Their status of "fossils" gives them a unique role 
to understand the origins of the Solar system, not only from the 
physical and dynamical point of view but also from the chemi- 
cal point of view (thanks to the knowledge of the compounds of 
the nucleus). When the comet gets closer to the Sun, the ices of 
the nucleus sublimate to form the coma (the comet atmosphere) 
where oxygen atoms are detected. Oxygen is an important ele- 
ment in the chemistry of the Solar system given its abundance 
and its presence in many molecules including H2O which con- 
stitutes 80% of the cometary ices. Oxygen atoms are produced 
by the photo-dissociation of molecules coming from the sub- 
limation of the cometary ices. Chemical reactions Eq.(l) to (6) 
involve possible parent molecules (Bhardwaj & Raghuram 2012; 
Festou & Feldman 1981). 



CO + hv 
CO + hv 



C + 0*CD) 
C + O'CS) 



(5) 
(6) 



Oxygen atoms have been detected in comets through the 
three forbidden lines observed in emission at 5577.339 A (the 
green line), 6300.304 A and 6363.776 A (the red lines) (Swings 
1962). These lines are coming from the deexcitation of the up- 
per state (2p 4 ) 1 S to the (2p 4 )'D 2 state for the 5577.339 A line 
and from the (2p 4 )'D2 state to the (2p 4 ) 3 Pi,2 states for the dou- 
blet red lines (see Fig. 1). The green line measurement is more 
difficult owing to its fainter intensity and the many C2 lines lo- 
cated around it. One of the first theoretical studies was carried 
out by Festou & Feldman (1981). They reviewed the production 
rate of O('S) and O('D) obtained in the laboratory from H2O, 
CO and CO2 photo-dissociations and measured the correspond- 
ing 'S/'D ratio. Recently, Bhardwaj & Raghuram (2012) made 
a new model for oxygen atom emissions and calculated the pro- 
duction rate of O('S) and O('D) for chemical reactions (1) to 
(6). The estimated 'S/'D ratios are significantly different from 
Festou & Feldman (1981) and are given in Table 1. 



H 2 + hv 
H 2 + hv 

C0 2 + hv 
C0 2 + hv 



H 2 + 0*( l D) 
H 2 + 0*CS) 

CO + 0\ l D) 
CO + 0*( l S) 



(1) 
(2) 

(3) 
(4) 



* Based on observations made with ESO Telescope at the La 
Silla Paranal Observatory under programs ID 268.C-5570, 270. C- 
5043, 073.C-0525, 274.C-5015, 075.C-0355, 080.C-0615, 280.C-5053, 
086.C-0958 and 087.C-0929. 



Up to now, little systematic researchwork has been done to 
study these lines at various heliocentric distances because the de- 
tection of the forbidden lines requires both high spectral and high 
spatial resolutions. Morrison et al. (1997) observed oxygen in 
comet C/1996 B2 (Hyakutake) with the 1-m Ritter Observatory 
telescope, Cochran & Cochran (2001) and Cochran (2008) anal- 
ysed [OI] lines in the spectra of 8 comets observed at the 
McDonald Observatory (see Bhardwaj & Raghuram (2012) for 
a complete review of these measurements). The present paper re- 
ports the results obtained for a homogeneous set of high quality 
spectra of 12 comets of various origins observed since 2003 with 
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(2p 4 ) 3 P„ rl , 2 



Fig. 1. Energy level diagram for [OI] lines. 



Parents 


Production rate (s ') 


Ratio 




0('S) 0('D) 


'S/'D 


'S/'D" 


H 2 


6.4 10~ 8 8.0 10-' 


0.080 


-0.1 


CO 


4.0 10~ 8 * 5.1 10~ 8 


0.784 


~1 


co 2 


7.2 10~ 7 1.2 10- 6 


0.600 


~1 



Table 1. The production rates and 'S/'D ratios obtained by Bhardwaj 
& Raghuram (2012) for photo-dissociation reactions at 1 AU. 
" The ratios obtained by Festou & Feldman (1981) are given in the last 
column for comparison. 

* This rate comes from Huebner & Carpenter (1979). 



UVES mounted on one of the 8-m telescope of the ESO VLT. 
First, we measured the intensities of the three forbidden oxygen 
lines. Two ratios were evaluated : the ratio between the two red 
lines, /630o/^6364> an d the ratio of the green line to the sum of the 
red lines, hsn I {Iaoa + hi6A) which we shall denote as G/R here- 
after. The purpose of the latter is to determine the main parent 
molecule of the oxygen atoms by comparing our results with the 
Bhardwaj & Raghuram (2012) effective excitation rates. We also 
observed how the G/R ratio depends on the heliocentric distance 
by analysing spectra of comets C/2001 Q4 (NEAT) and C/2009 
PI (Garradd) at small and large distances from the Sun. 
Finally, we measured the Full Width at Half Maximum (FWHM) 
of the three lines. There is a long-standing debate on the FWHM 
measurement because the first observations made by Cochran 
(2008) are in contradiction with the theory. The intrinsic width 
of the green line is wider than the red ones while the theory pre- 
dicts the opposite (Festou 1981). 

2. Observations 

Our analysis is based on 12 comets listed in Table 2, consider- 
ing the two components of 73P/Schwassmann-Wachmann 3 (B 
and C) as two comets. This sample is characterized by a large 
diversity. These comets have various dynamical origins (exter- 
nal, new, Jupiter family, Halley type) and are observed at dif- 
ferent heliocentric distances (from 0.68 AU to 3.73 AU). Some 
are split comets like 73P/SW3 or did a close approach to Earth. 
Others have been observed at large heliocentric distances like 
Q4 (NEAT) and PI (Garradd). The observing material is made 
of a selection of 48 high signal-to-noise spectra obtained during 
the last ten years (Manfroid et al. (2009) ; Jehin et al. (2009) 
and references therein) with the cross-dispersed echelle spectro- 
graph UVES (Dekker et al. 2000). A dichroic splits the light 
beam into two arms (3000 to 5000 A and 4200 to 1 1000 A). In 
most observations, a slit of 0.45" x 11" was used. With such a 
slit, the resolving power is R- A/AA = 1 10000 in the red arm. 



All spectra were recorded with nearly the same instrumental set- 
ting. Such a high resolution is needed to isolate the [OI] emis- 
sion lines from the telluric lines and other cometary lines (C2, 
NH2) (see Fig. 2). The slit was always centered on the cometary 
nucleus. For a comet located at 1 AU from the Earth, the spa- 
tial area covered by the slit is approximatively 320 km by 8000 
km. The size of the slit for each individual spectrum is given in 
Table 3. 

3. Data reduction 

The 2D spectra were reduced with the UVES pipeline reduction 
program. For data obtained until 2008, the UVES pipeline (ver- 
sion 2.8.0) was used within the ESO-MIDAS environment to ex- 
tract, for each order separately, the 2D spectra, bias-subtracted, 
flat-fielded and wavelength calibrated. ESO-MIDAS is a spe- 
cial software containing packages to reduce ESO data 1 . For a 
given setting, the orders were then merged using a weighting 
scheme to correct for the blaze function, computed from high 
S/N master flat-fields reduced in exactly the same way. This pro- 
cedure leads to a good order merging. The post-2008 observa- 
tions were reduced with the gasgano CPL 2 interface (version 
2.4.3) of the re-furbished UVES pipeline that directly provides 
accurately merged 2D spectra. ID spectra were extracted by av- 
eraging the 2D ones, with cosmic rejection. Standard stars were 
similarly reduced and combined to derive instrumental response 
functions used to correct the ID spectra. The wavelength cal- 
ibration was computed using Thorium Argon (Th-Ar) spectra. 
The calibration revealed small shifts of the lines due to the fact 
that the Th-Ar spectra were usually taken in the day and not 
right after the comet observations. This shift was removed using 
9 cometary lines of NH2 in the vicinity of the red [OI] lines and 
for which laboratory wavelengths are well known. 
Before doing any measurement, the telluric absorption lines 
were removed and the solar continuum contribution subtracted 
using the BAAS and Kurucz spectra. These two spectra are so- 
lar spectra, with atmospheric lines for the BAAS 3 spectrum and 
without for the Kurucz 4 one (2005). Around the red doublet line 
at 6300 A, there are indeed telluric absorption lines mostly due 
to O2 molecules in the Earth's atmosphere which could lead to an 
underestimate of the forbidden oxygen line intensity (see Fig. 3). 

The Solar continuum has three contributions. First and fore- 
most is the one due to the reflection of sunlight on the dust parti- 
cles in the coma. Another one appears when the observations are 
made in or close to the twilight. A third contribution is, in a few 
cases, the background radiation by the Moon. To remove these 
contaminations, we Doppler shifted the BAAS spectrum to the 
proper values and scaled the intensities until the solar features 
are completely removed. This provides the final spectrum such 
as the one presented in Fig. 4 for comet 8P/Tuttle, free of any 
solar and telluric lines. 

4. Data analysis 

After the reduction and the correction of the data, we measured 
both the intensity and the FWHM of the three forbidden oxy- 
gen lines by fitting a gaussian profile using the IRAF 5 software. 

1 http://www.eso.org/sci/software/esomidas/ 

2 http://www.eso.org/sci/software/gasgano/ 

3 http://bass2000.obspm.fr/solar_spect.php 

4 http://kurucz.havard.edu/sun/irradiance2005/irradthu.dat 

5 IRAF is a tool for the reduction and the analysis of astronomical 
data (http://iraf.noao.edu). 
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Comet 


Short name comet 


Tp 


e 


a 


q 


i 


P 


Tj 


Type (L) 


C/2002 VI (NEAT) 


VI (NEAT) 


18-02-2003 


0.99990 


1010 


0.10 


82 


32100 


0.06 


EXT 


C/2002 X5 (Kudo-Fujikawa) 


X5 (Kudo-Fujikawa) 


29-01-2003 


0.99984 


1175 


0.19 


94 


40300 


-0.03 


EXT 


C/2002 Yl (Juels-Holvorcem) 


Yl (Juels-Holvorcem) 


13-04-2003 


0.99715 


250.6 


0.71 


104 


3967 


-0.23 


EXT 


C/2001 Q4 (NEAT) 


Q4 (NEAT) 


16-05-2004 


1.00069 




0.96 


100 






NEW 


C/2002 T7 (LINEAR) 


T7 (LINEAR) 


23-04-2004 


1.00048 




0.61 


161 


_ 




NEW 


C/2003 K4 (LINEAR) 


K4 (LINEAR) 


14-10-2004 


1.00030 


_ 


1.02 


134 


_ 


_ 


NEW 


9P/Tempel 1 


Tempel 1 


05-07-2005 


0.51749 


3.1 


1.51 


11 


5.5 


2.97 


JF 


73P-C/Schwassmann-Wachmann 3 


SWc3 


07-06-2006 


0.69338 


3.1 


0.94 


11 


5.4 


2.78 


JF 


73P-B/Schwassmann-Wachmann 3 


SWb3 


08-06-2006 


0.69350 


3.1 


0.94 


11 


5.4 


2.78 


JF 


8P/Tuttle 


Tuttle 


27-01-2008 


0.8198 


5.7 


1.03 


55 


13.6 


1.60 


HF 


103P/Hartley 2 


Hartley 2 


28-10-2010 


0.6950 


3.5 


1.06 


14 


6.5 


2.64 


JF 


C/2009 PI (Garradd) 


PI (Garradd) 


23-12-2011 


1.0011 




1.55 


106 






NEW 



Table 2. Orbital characteristics of the 12 comets and their classification. T p is the epoch of the perihelion (dd-mm-yyyy), e the eccentricity, a is related to the 
semi-major axis of the orbit (AU), q is the perihelion distance (AU), i represents the inclination on the ecliptic (degrees), P is the period (years), Tj is the Tisserand 
parameter relative to Jupiter. The last data correspond to the classification of the comets according to Levison (1996) : HF and JF mean "Halley Family" and "Jupiter 
family" and gather comets with short period (< 200 years) ; EXT and NEW are respectively for external comets (a<10000 AU) and new comets (a>10()00 AU) 
which directly come from the Oort Cloud. 
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3 329.344 


1.20 


14.81 


1.51 


-28.23 
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12.00 


9P/Tempel 1 


3 553.955 


1.51 


-0.21 


0.89 


8.95 


7200 


0.44 x 


12.00 


9P/Tempel 1 


3 554.954 


1.51 


-0.15 


0.89 


9.07 


7200 


0.44 x 


12.00 


9P/Tempel 1 


3 555.955 


1.51 


-0.04 


0.90 


9.19 


7200 


0.44 x 


12.00 


9P/Tempel 1 


3 557.007 


1.51 


0.09 


0.90 


9.48 


9600 


0.44 x 


12.00 


9P/Tempel 1 


3 557.955 


1.51 


0.20 


0.91 


9.44 


7500 


0.44 x 


12.00 


9P/Tempel 1 


3 558.952 


1.51 


0.31 


0.91 


9.55 


7500 


0.44 x 


12.00 


9P/Tempel 1 


3 559.954 


1.51 


0.43 


0.92 


9.68 


7500 


0.44 x 


12.00 


9P/Tempel 1 


3 560.952 


1.51 


0.55 


0.93 


9.80 


7800 


0.44 x 


12.00 


9P/Tempel 1 


3 561.953 


1.51 


0.66 


0.93 


9.91 


7200 


0.44 x 


12.00 


9P/Tempel 1 


3 562.956 


1.51 


0.78 


0.94 


10.04 


7200 


0.44 x 


12.00 


73P-C/SW 3 


3 882.367 


0.95 


-4.17 


0.15 


12.31 


4800 


0.60 x 


12.00 


73P-B/SW 3 


3 898.369 


0.94 


1.79 


0.25 


13.10 


4800 


0.60 x 


12.00 


8P/Tuttle 


4 481.021 


1.04 


-4.29 


0.36 


21.64 


3600 


0.44 x 


10.00 


8P/Tuttle 


4 493.018 


1.03 


0.40 


0.52 


24.72 


3900 


0.44 x 


10.00 


8P/Tuttle 


4 500.017 


1.03 


3.16 


0.62 


24.16 


3900 


0.44 x 


10.00 


103P/Hartley 2 


5 505.288 


1.06 


2.53 


0.16 


7.08 


2900 


0.44 x 


12.00 


103P/Hartley 2 


5 505.328 


1.06 


2.55 


0.16 


7.19 


3200 


0.44 x 


12.00 


103P/Hartley 2 


5 510.287 


1.07 


4.05 


0.18 


7.96 


2900 


0.44 x 


12.00 


103P/Hartley 2 


5 510.328 


1.07 


4.07 


0.18 


8.07 


3200 


0.44 x 


12.00 


103P/Hartley 2 


5 511.363 


1.08 


4.37 


0.19 


8.27 


900 


0.44 x 


12.00 


C/2009 PI (Garradd) 


5 692.383 


3.25 


-16.91 


3.50 


-44.66 


3600 


0.44 x 


12.00 


C/2009 PI (Garradd) 


5 727.322 


2.90 


-16.89 


2.57 


-46.38 


3600 


0.44 x 


12.00 


C/2009 PI (Garradd) 


5 767.278 


2.52 


-16.46 


1.64 


-29.26 


1800 


0.44 x 


12.00 


C/2009 PI (Garradd) 


5 813.991 


2.09 


-14.82 


1.47 


14.79 


4800 


0.44 x 


12.00 


C/2009 PI (Garradd) 


5 814.974 


2.08 


-14.77 


1.48 


15.31 


4800 


0.44 x 


12.00 


C/2009 PI (Garradd) 


5 815.982 


2.07 


-14.71 


1.49 


15.84 


4800 


0.44 x 


12.00 



Table 3. Individual spectra. MJD = JD - 2 450 000.5 the Modified Julian Day, r is the heliocentric distance, A the geocentric distance, r and A are respectively 
the heliocentric and geocentric velocities. Exptime corresponds to the exposure time in seconds. Slit gives the size of the entrance slit of the spectrograph in arc 
seconds. 
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Fig. 3. The telluric absorption line correction in the vicinity of red [OI] 
lines. The black spectrum has been corrected (thick line), the red one is 
uncorrected. 
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Fig. 4. The solar continuum subtraction in the vicinity of the red [OI] 
lines. The solar continuum has been removed in the lower spectrum 
(black) while it is present in the upper one (in red). It corresponds to the 
final spectrum on which the measurements of the [OI] will be made. 



FWHM,,! 



FWHM . 
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Fig. 2. The 6300.304 A, 6363.776 A and 5577.339 A oxygen lines (tel- 
luric (©) and cometary) in a 3600s spectrum of comet 8P/Tuttle ob- 
tained in March 2008. Thanks to the high resolution of the UVES spec- 
trograph, the cometary and telluric lines are well separated when the 
Doppler shift is larger than 15 km s -1 . 



The observed width is the convolution of the instrumental profile 
with the natural width : 



where the instrumental width corresponds to the width of the 
Th-Ar lines. 



These measurements can only be performed with best accuracy 
when the cometary oxygen lines are well separated from the tel- 
luric [OI] lines, i.e, when the geocentric velocities of the comet 
exceeds 15 km s _1 . Using the deblend function of the splot pack- 
age in IRAF, we could measure the [OI] lines when the Doppler 
shift is as low as 7 km s . 
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Comet 


MJD 


r(AU) 


^6300 
6364 


'5577 
^6300+^6364 




FWHM mlrmslc 














5577.339 A 


6300.304 A 


6363.776 A 


C/2002 VI (NEAT) 


2 647.037 


1.22 


3.10 


0.06 


2.08 


1.57 


1.40 


C/2002 VI (NEAT) 


2 647.062 


1.22 


3.02 


0.06 


1.95 


1.42 


1.44 


C/2002 VI (NEAT) 


2 649.031 


1.18 


3.05 


0.07 


2.23 


1.45 


1.42 


C/2002 VI (NEAT) 


2 649.056 


1.18 


3.04 


0.07 


2.11 


1.40 


1.38 


C/2002 VI (NEAT) 


2 719.985 


1.01 


3.05 


0.10 


2.76 


1.54 


1.49 


C/2002 X5 (Kudo-Fujikawa) 


2 705.017 


1.06 


2.97 


0.10 


2.64 


1.46 


1.50 


C/2002 X5 (Kudo-Fujikawa) 


2 705.039 


1.07 


2.98 


0.10 


2.61 


1.54 


1.60 


C/2002 X5 (Kudo-Fujikawa) 


2 705.060 


1.07 


3.01 


0.10 


2.64 


1.50 


1.54 


C/2002 Yl (Juels-Holvorcem) 


2 788.395 


1.14 


2.93 


0.09 


2.88 


1.52 


1.47 


C/2002 Yl (Juels-Holvorcem) 


2 788.416 


1.14 


3.01 


0.09 


2.75 


1.54 


1.49 


C/2002 Yl (Juels-Holvorcem) 


2 789.394 


1.16 


2.98 


0.11 


3.09 


1.59 


1.57 


C/2002 Yl (Juels-Holvorcem) 


2 789.415 


1.16 


3.01 


0.10 


3.03 


1.61 


1.59 


C/2001 Q4 (NEAT) 


2 883.293 


3.73 


2.95 


0.33 


2.51 


2.39 


2.39 


C/2001 Q4 (NEAT) 


2 883.349 


3.73 


3.22 


0.33 


2.31 


2.63 


2.69 


C/2001 Q4 (NEAT) 


2 889.236 


3.67 


3.17 






2.49 


2.52 


C/2001 Q4 (NEAT) 


2 889.320 


3.66 


3.51 


0.30 


2.55 


2.75 


2.36 


C/2002 T7 (LINEAR) 


3 131.421 


0.68 


3.04 


0.14 


3.12 


1.82 


1.80 


C/2002 T7 (LINEAR) 


3 151.976 


0.94 


3.51 


0.11 


2.84 


1.70 


1.71 


C/2002 T7 (LINEAR) 


3 152.036 


0.94 


3.39 


0.09 


2.81 


1.66 


1.72 


C/2003 K4 (LINEAR) 


3 131.342 


2.61 


3.08 


0.10 


2.51 


1.81 


1.66 


C/2003 K4 (LINEAR) 


3 132.343 


2.59 


3.13 


0.08 


2.76 


2.12 


1.94 


C/2003 K4 (LINEAR) 


3 329.344 


1.20 


3.20 


0.07 


2.38 


1.52 


1.52 


9P/Tempel 1 


3 553.955 


1.51 


3.06 


0.05 


2.91 


1.61 


1.43 


9P/Tempel 1 


3 554.954 


1.51 


3.00 


0.04 


2.14 


1.59 


1.48 


9P/Tempel 1 


3 555.955 


1.51 


3.00 






1.56 


1.46 


9P/Tempel 1 


3 557.007 


1.51 


3.05 


0.05 


2.03 


1.55 


1.44 


9P/Tempel 1 


3 557.955 


1.51 


3.23 


0.03 


1.75 


1.54 


1.47 


9P/Tempel 1 


3 558.952 


1.51 


3.23 


0.04 


1.79 


1.58 


1.48 


9P/Tempel 1 


3 559.954 


1.51 


3.21 


0.05 


2.57 


1.84 


1.60 


9P/Tempel 1 


3 560.952 


1.51 


3.17 


0.06 


2.77 


1.72 


1.52 


9P/Tempel 1 


3 561.953 


1.51 


3.18 


0.04 


1.92 


1.63 


1.48 


9P/Tempel 1 


3 562.956 


1.51 


3.25 


0.03 


1.67 


1.61 


1.48 


73P-C/Schwassmann-Wachmann 3 


3 882.367 


0.95 


3.24 


0.11 


2.12 


1.31 


1.30 


73P-B/Schwassmann-Wachmann 3 


3 898.369 


0.94 


3.03 


0.11 


2.46 


1.54 


1.48 


8P/Tuttle 


4 481.021 


1.04 


3.04 


0.05 


2.16 


1.55 


1.36 


8P/Tuttle 


4 493.018 


1.03 


3.03 


0.04 


2.21 


1.55 


1.41 


8P/Tuttle 


4 500.017 


1.03 


3.03 


0.05 


2.28 


1.47 


1.33 


103P/Hartley 2 


5 505.288 


1.06 


3.12 


0.09 


2.30 


1.35 


1.29 


103P/Hartley 2 


5 505.328 


1.06 


3.22 


0.10 


2.31 


1.26 


1.25 


103P/Hartley 2 


5 510.287 


1.07 


3.24 


0.07 


2.31 


1.26 


1.06 


103P/Hartley 2 


5 510.328 


1.07 


3.17 


0.08 


2.24 


1.27 


1.17 


103P/Hartley 2 


5 511.363 


1.08 


3.17 


0.09 


2.19 


1.23 


1.18 


C/2009 PI (Gairadd) 


5 692.383 


3.25 


2.92 


0.21 


2.28 


1.67 


1.87 


C/2009 PI (Garradd) 


5 727.322 


2.90 


3.38 


0.14 


2.54 


1.27 


1.49 


C/2009 PI (Gairadd) 


5 767.278 


2.52 


3.12 


0.10 


2.19 


1.30 


1.31 


C/2009 PI (Gairadd) 


5 813.991 


2.09 


3.01 


0.07 


2.18 


1.25 


1.33 


C/2009 PI (Garradd) 


5 814.974 


2.08 


3.07 


0.07 


2.16 


1.33 


1.33 


C/2009 PI (Garradd) 


5 815.982 


2.07 


3.04 


0.06 


2.21 


1.57 


1.58 



Table 4. Measured line ratios and velocity widths (km s ) for the three forbidden oxygen lines in the spectra of each comet. For 
the third spectrum of comets Q4 (NEAT) and Tempel 1, the green line couldn't be used because of a contamination by a cosmic ray. 



5. Results 

5.1. Intensity ratios 

The two measured intensity ratios, hm)lhiM and 
G/R = hsnlilaoo + ^6364 X are displayed for each spec- 
trum in Fig. 5 and Fig. 6 respectively. When the collisional 
quenching is neglected, the intensity of a line is given by (Festou 
&Feldman 1981): 

/ = 1(T 6 r p l a /3 N photons cm" 2 s" 1 (8) 

where t p is the photo-dissociative lifetime of the parent, a is the 
yield of photo-dissociation, f3 corresponds to the branching ra- 
tio for the transition and /V is the column density of the parent. 
Table 4 lists the results of the ratios for the 12 comets. The red 
doublet ratio is remarkably similar for all comets. The average 
value over the whole sample of spectra is 3.12 + 0.02. The error 
corresponds to the average sample standard deviation (cr/ VAO. 



For a given comet, the dispersion of the ratio measurements from 
one spectrum to another is about 2-3%. This result is in very 
good agreement with the theoretical values of the branching ra- 
tio given by quantum mechanics, 3.096 by Galavis et al. (1997) 
and 2.997 by Storey & Zeippen (2000). Since both red lines are 
transitions from the same level to the ground state, t p , a and N 
in Eq. (8) are equal and the intensity ratio /630o/^6364 is indeed 
equivalent to the branching ratio fiscal PaaA- The ratio is also in 
agreement with the value of 3.09 + 0.12 obtained by Cochran 
(2008) from a sample of 8 comets. 

The G/R ratio has an average value of 0. 1 1 + 0.07 for the 
complete sample (see Fig. 6). The error corresponds to the sam- 
ple standard deviation (cr) and will be given for the other aver- 
age measurements. The dispersion of the ratio measurements for 
all spectra of a given comet is about 10 - 15%. This result is 
in agreement with the value of 0.09 + 0.01 found by Cochran 
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(2008). This leads to the conclusion that H2O is the main par- 
ent molecule producing oxygen atoms according to Bhardwaj & 
Raghuram (2012) values (Table 1). 
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Fig. 5. The red doublet ratio for each comet. The spectra of a given 
comet have the same symbols. The average value of the sample 
(3.12 ± 0.11) is shown by the solid line. The dashed lines show the 
standard deviation (cr). 
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Fig. 6. Same figure as Fig. 5 for the G/R ratio. The mean value is 
0.1 1 ± 0.07. Note the data points for comet Q4 (NEAT) at 3.7 AU. 



5.2. Line widths 

Table 4 lists the intrinsic line velocity widths (FWHM (v)). The 
FWHM (v) is obtained from the FWHM intrinsic (T) given in the 
Eq. 7 thanks to the relation : 



FWHM (v) 



FWHM;, 



(9) 



i[orj 



where A[oi] corresponds to the wavelength of the considered 
oxygen line. The dispersion of the FWHM (v) for a given comet 
is about 5% for the red lines and 10% for the green line which 
is fainter. The two red lines widths are equal within the errors 
(1.61 ± 0.34 km s" 1 for the 6300 A and 1.56 ± 0.54 km s" 1 
for the 6364 A line) which is not a surprise since both lines 
are transitions from the O('D) state to the ground state (Fig. 7). 
This result is also consistent with the values of Cochran (2008) 
(1.22 ± 0.36 km s" 1 ). Figs. 8, 9 and 10 present the width of 
the three [OI] cometary lines, the three [OI] telluric lines and 
the width of some representative cometary lines in the vicin- 
ity. Those lines correspond to NH2 lines in the red region and 
to C2 lines in the green region. The widths of the red telluric 
lines are slightly smaller than the [OI] cometary lines (1.37 ± 
0.35 km s _1 for the 6300.304 A line, 1.36 ± 0.39 km s~' for 
the 6363.776 A line) but their values are in quite good agree- 
ment with the cometary lines. Fig. 1 1 shows clearly that the [OI] 
cometary green lines are wider than the two red lines. Our mean 
value for the [OI] cometary green line is 2.44 + 0.28 km s -1 (to 
be compared to 0.94 ± 0. 16 km s _1 for C 2 lines and 0.90 + 0.37 
km s _1 for the green telluric line). This peculiarity was al- 
ready noticed by Cochran (2008) who found a mean velocity 
of 2.49 ± 0.36 km s _1 in good agreement with our value. This 
width could be explained if the excess energy for the formation 
of the O('S) state is larger than for the O('D) state. Contrary to 
what is observed, theoretical models using Ly-a photons as exci- 
tation source give an excess velocity of 1.6 km s for the O('S) 
state and a large value of 1.8 km s for O('D) in the case of 
water photo-dissociation (Festou 1981). 



6. Discussion 

6.1. The G/R ratio at large heliocentric distance 

In Fig. 6, we noticed that C/2001 Q4 (NEAT) has a G/R ratio 
of about 0.3 clearly such higher than other comets with values 
around 0.1. Q4 (NEAT) was at 3.7 AU from the Sun and this 
peculiarity led us to look more carefully on comets observed at 
large heliocentric distances (Fig. 12) (Decock et al. 201 1). 

If we only consider the data taken at r < 2 AU, the G/R ratio 
average value is equal to 0.09 ± 0.02 which is in good agreement 
with the ratio obtained by Bhardwaj & Raghuram (2012) and 
Festou & Feldman (1981) for H2O as the parent molecule (see 
Table 1). The large values of G/R in Q4 (NEAT) at 3.7 AU could 
be explained by the increasing contribution, at large heliocentric 
distances, of other parent molecules producing oxygen atoms. 
It is well known that the sublimation of H2O significantly de- 
creases beyond 3 AU while the sublimation of other ices like CO 
and/or CO2 dominates comet activity (i.e., Crovisier & Encrenaz 
2000). In order to investigate this hypothesis, we observed comet 
C/2009 PI (Garradd) at four different heliocentric distances from 
3.25 AU to 2.07 AU. Fig. 13 presents the G/R intensity ratios as 
a function of r for all spectra. A relation between the heliocen- 
tric distance and the G/R intensity ratio can be seen : the latter is 
getting rapidly larger when the comet is getting farther from the 
Sun. We added in this plot the Cochran & Cochran (2001) and 
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Fig. 7. The 6300.304 A line width against the 6363.776 A one. The 
different spectra for a given comet are denoted with the same symbol. 
Note the larger width for Q4 (NEAT). 



Cochran (2008) average values coming from their sample of 8 
comets. All these values are clumped around 1 AU with a ratio 
of ~0. 1 . Q4 (NEAT) is one of their 8 comets and its G/R ratio is 
0.09 + 0.02 at r = 0.98 AU. A few other measurements found 
in the literature for comets observed at heliocentric distances be- 
yond 2 AU were also added to the graph. Furusho et al. (2006) 
used the Subaru telescope to analyse the forbidden oxygen lines 
in comet 116P/Wild 4 observed at r — 2.4 AU. Capria et al. 
(2010) measured the intensity of these lines during and after the 
outburst of comet 17P/Holmes which occurred on October 24, 
2007 when the comet was located at about r = 2.5 AU. McKay 
et al. (2012b) studied the [OI] lines of the two comets C/2006 
W3 Christensen and C/2007 Q3 Siding Spring, respectively at 
3.13 and 2.96 AU, with the ARCES echelle spectrometer of 3.5- 
m telescope at Apache Point Observatory. While the quality of 
the data are not as good because the lines are faint and sometimes 
heavily blended with the telluric line, all these measurements 
show a relatively high value of the G/R ratio when r > 3 AU. 
This confirms the hypothesis that oxygen is also coming from 
other molecules. CO and/or CO2 are obvious candidates as they 
produce large values of the G/R ratio (see Table 1). 

Using HST/STIS 6 observations, Feldman et al. (2004) esti- 
mated at 4% the CO abundance relative to water in comet Q4 
(NEAT) located at 1 AU from the Sun. The study made by Biver 
et al. (2012) on comet PI (Garradd) provided an average value of 
5% for the CO abundance when the comet was at 1.9 AU at pre- 
perihelion phase (October 201 1). At large heliocentric distances, 
CO/H2O had to be even higher. So the CO molecule could be a 
candidate to explain such a large G/R ratio. However, Bhardwaj 
& Raghuram (2012) showed that after the photo-dissociation of 
H2O, the next main source for the green line emission is CO2 
with 10-40% abundance relative to water. They compared their 
model with the observational values obtained for C/1996 B2 



Fig. 8. The width velocity of the 6300.304 A [OI] cometary line (■), the 
[OI] telluric line (a) and nearby NH 2 cometary lines (•) (km s~') for 
one spectrum per comet. The typical errors correspond to the standard 
deviation (a) obtained considering all the spectra of comet 8P/Tuttle. 
We have also checked this error by measuring randomly <x for a few 
spectra of other comets. For Q4 (NEAT) and K4 (LINEAR), the NH 2 
cometary lines are missing because no lines were detected in this range. 
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6 HST/STIS is the acronym for Hubble Space Telescope / Space 
Telescope Imaging Spectrograph (http://www.stsci.edu/hst/stis). 



Fig. 9. Same figure as Fig. 8 for the 6363.776 A line. 



(Hyakutake) by Morrison et al. (1997) and Cochran (2008) when 
the comet was around 1 AU. Hyakutake is also rich in CO and 
its abundance was evaluated at 22% at 1 AU (Biver et al. 1999). 
Assuming a CO2 abundance of 1% and a O('S) yield of 0.2%, 
Bhardwaj & Raghuram (2012) concluded that the production 
rates of the O('S) are similar for CO and C0 2 (-10-25%) with 
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Fig. 10. Same figure as Fig. 8 for the 5577.339 A line and with C2 



cometary lines. The larger width, by about 
cometary green line is obvious. 
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Fig. 11. The 5577.339 A line width against the average width of the 
red lines. The different spectra of a given comet are denoted with the 
same symbol. Note that the green line is always larger than the red lines 
except for Q4 (NEAT). 



such abundances. Therefore, since Q4 (NEAT) and PI (Garradd) 
have less CO than Hyakutake, CO2 molecules could rather be the 
main contributor to the formation of the green line emission for 
these two comets observed at large heliocentric distance. CO2 
measurements are unfortunately very rare because they are dif- 
ficult to obtain to confirm this hypothesis. But if this is the case, 
G/R ratios estimated might be used as a proxy for the CO2 rela- 
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Fig. 12. These two hours spectra of comet Q4 (NEAT) obtained with 
the ESO VLT are the most distant detection of [OI] lines in a comet (at 
3.7 AU). The lines are well detected and well separated from the telluric 
line. The width of the green and the red lines are the same. 



tive abundance but such a relation should be carefully calibrated 
on a sample of comets with known CO2, as also suggested by 
McKay et al. (2012a). 
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Fig. 13. The G/R intensity ratio as a function of the heliocentric dis- 
tance (AU). The same symbol is used for the different spectra of each 
comet. The most distant comets have a higher ratio. Other values from 
the literature are also plotted. The solid line represents the average value 
only for our sample of spectra and the standard deviation (cr) range is 
denoted with horizontal dashed lines. If we only consider the data at 
r < 2 AU, the average value of G/R ratio is 0.09 ± 0.02. The vertical 
dashed line at 3 AU represents the distance out to which the sublimation 
of water is strongly decreasing (Crovisier & Encrenaz 2000). 



Using Eq. 8 and only considering the main parent species of 
oxygen atoms (i.e H2O, CO2 and CO), we can write the G/R 
ratio as : 



The W b are independent of the heliocentric distance because 
all of them depend on the same heliocentric distance. This is 
also the case for the values given in Table 1. Considering the 
same values for and 0« reen ! we also computed the CO2/H2O 
abundance from the G/R ratios of comets Q4 (NEAT) and 
PI (Garradd) measured by Cochran (2008) and McKay et al. 
(2012a), respectively. From our data, we find an exponential 
relation between the CO2/H2O abundance and the heliocentric 
distance of the type C0 2 /H 2 = A e br with A = 1 10 4 and 
b = 2.41. These results are shown in Fig. 14. We noticed that 
CO2/H2O is as high as 75% in Q4 (NEAT) and that C0 2 starts 
to contribute to the G/R ratio at 2.5 AU in the comets. 
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preen W%% /3 S '"" N H ,0 + WgT P"'"" N CO , + W S ™ /}*'"" N CO 

G/R = — -r = — : : h : : • (1°) 

jred W red g pet ff^ + fl^rf fired ff^ + W * fired ^ 

where W£ = r^ 1 corresponds to the effective produc- 
tion rates of O('S) and O('D) states 7 . However, if we consider 
that oxygen atoms only come from H2O and CO molecules, 
we found that G/R ratio does not change strongly : it varies 
from 0.08 to 0.1 for a CO abundance between 10% and 80%. 
Therefore, we assume that the [OI] atoms are only produced by 
the H2O and CO2 photo-dissociations. Hence, Eq. 10 can be sim- 
plified as : 



Fig. 14. Evolution of the CO2/H2O abundance as a function of the 
heliocentric distance (AU). Assuming that a bigger G/R ratio comes 
from a more important contribution of C0 2 for large heliocentric dis- 
tance comets, measuring this ratio can give us the CO2 abundance from 
ground observations. The curve corresponds to an exponential fit of our 
data (■ and •). 

Therefore, the forbidden oxygen lines measurement could 
become a new way to determine the CO2 abundance in comets at 
different heliocentric distances from ground observations while 
the direct measurement of CO2 molecules is only possible from 
space. 



jgreen Wl™ P green Nh 2 + ^rT P*™" N C 2 

G/R = — t = 1 z : ^7 ; (11) 

ped W red o yed ^ + W ^ pred ^ 

Therefore, the CO2/H2O abundance can be computed by : 

(12) 



Nc0i _ (G/R) W r H % -^ ree " WgJ 

N Hl O ~ fiSreen yjZ™™ _ (Q /R) W gd 



The production rates W? are given in Tables 1 and 2 of Bhardwaj 



& Raghuram (2012) (i.e W™ = 8 10"' s" 



W g ™ n = 7.2 10- 



H 2 

, W$. = 1.2 10- 6 s" 1 , W s Z e " = 1 



6.4 10- s s" 



WJ!i = 5.1 10" 8 s -1 ) and p s " e " = 0.91 (Slanger et al. 2006). 



6.2. The [OI] lines widths 

One explanation for the large width of the [OI] green line could 
be that exciting photons producing oxygen atoms in the ! S 
state have more energy than the Lyman a photons and/or come 
from an other parent molecule. Indeed, Fig. 7 of Bhardwaj & 
Raghuram (2012) shows that the wavelength region for the ma- 
jor production of O('S) coming from CO2 seems to be the 955- 
1 165 A band. New theoretical studies could be necessary to eval- 
uate the excess velocity for the O('S) state coming from the 
photo-dissociations of both H2O and CO2. 

In Fig. 1 1, we also notice that the [OI] red lines of the comet 
Q4 (NEAT) are wider than in other comets while the width 
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G/R 


A/^.rt /A/>, 
iyj CCh I^HiO 


5 505.288 


1.06 


0.09 


0.028 


5 505.328 


1.06 


0.10 


0.040 


5 510.287 


1.07 


0.07 


0.001 


5 510.328 


1.07 


0.08 


0.015 


5 511.363 


1.08 


0.09 


0.030 



Table 5. The C0 2 /H 2 abundances estimated from Eq. 12 for the comet 
103P/Hartley 2. 



of the green line is normal with respect to other comets. This 
peculiarity seems real as both red lines show the same width 
and they are 2-cr away from the average value of the sample 
(1.58 + 0.30 km s~'). Moreover, the telluric lines in Q4 (NEAT) 
have the usual width (0.99 ± 0.10 km s _1 for the 6300 A line, 
0.87 ± 0.04 km s _I for the 6364 A line and 0.80 + 0.26 km s~' 
for the green line) which discards a problem with the spectra 
or the analysis. The width of the red lines and the green line is 
similar in this unique case (2.57 + 0.16 km s _1 for the 6300 A 
line, 2.49 ± 0.15 km s" 1 for the 6364 A and 2.46 ± 0.13 km s" 1 
for the 5577 A line) and could give us clues about the process 
at play. As previously discussed, at large heliocentric distances, 
the [OI] lines could mainly come from the CO2 molecules which 
preferentially photo-dissociate in the 1 S state while at low helio- 
centric distances, they essentially come from H2O which pref- 
erentially dissociates in the 'D state (Table 1). This means that 
at large distances, the red lines are mainly produced through the 
'S-'D channel which also produces the green line, while at low 
distances as H2O dominates they are mostly produced directly 
from the : D state. We thus expect the widths of the red and 
green lines to be similar at high distances because the O('S) and 
O('D) atoms are mostly produced from a same molecule, while 
they can differ at low distances, as observed. CO2 is the best 
candidate and the larger width could be explained by the main 
excitation source of CO2, the 955-1165 A band, which is more 
energetic than Ly-cr photons, as shown in Fig. 1 1 of Bhardwaj 
& Raghuram (2012). This particularity is not seen for comet 
PI (Garradd) maybe because the spectrum is not taken at suffi- 
ciently large heliocentric distance and CO2 abundance was lower 
than in Q4 (NEAT). Anyway, this hypothesis needs to be con- 
firmed by observing other comets at large heliocentric distances 
(> 3.5 AU). 

6.3. 103P/Hartley2 

The Jupiter Family comet 103P/Hartley 2 has been found rich in 
C0 2 (~ 20%) from EPOXI measurements (AHearn et al. 201 1) 
and poor in CO (Weaver et al. 201 1). Thanks to the Eq. 12, we 
could evaluate the CO2 abundances from our G/R ratios mea- 
sured at ~ 1 AU and compare them with the values obtained by 
EPOXI. The results are given in Table 5. The mean G/R ratio of 
0.09 ± 0.01 is similar to those of other comets observed below 
2.5 AU and give a relative abundance of CO2 of only ~ 3% while 
a G/R ratio of ~ 0.2 would be needed to reach a CO2 abundance 
as high as ~ 20%. 

We do not confirm a high value for 103P/Hartley 2 as 
claimed by (McKay et al. 2012) or a large dispersion in the abun- 
dances. Their values might be higher because they used a W^™" 

of Bhardwaj & Raghuram (2012) equal to 2.6 10~ 8 s _1 while the 
value provided in Table 1 of this paper is 6.4 10~ 8 s _1 . 

The discrepancy with the EPOXI CO2 abundance is inter- 
esting. Below 2.5 AU, G/R values are distributed from 0.05 for 
comets like 9P/Tempel 1 and 8P/Tuttle and can go up to 0. 12 for 



T7 (LINEAR). The spread from comet to comet is higher than 
the errors estimated from the dispersion from the spectra of a 
same comet which could be explained by an intrinsic variation 
of G/R from comet to comet due to the content of CO2. The G/R 
range, of the order of 0. 1 (implying a CO2 variation of ~ 20%) 
is of the same order as the range of CO2 abundance in comets of 
various origins measured below 2.5 AU (A'Hearn 2012). 

The fact that some comets, like 9P/Tempel 1 and 8P/Tuttle, 
have values lower than the pure water case (giving a 
G/R min = 0.08 from Bhardwaj & Raghuram (2012)) and that 
the 103P/Harley 2 G/R ratio is too low compared to EPOXI 
CO2 abundances, could indicate a problem with the models. 
The pure water ratio should be 0.05 or lower based on the 
comets with the smallest values (9P/Tempel 1 has a CO2 abun- 
dance of ~ 7% (Feaga et al. (2007a), Feaga et al. (2007b)). 
Bhardwaj & Raghuram (2012) tried values from 0% to 

1% and finally chose 1%, but this value is uncertain. Without 
C0 2 , G/R min = j3 green W s ™ e n I W™* . If we use 0.5% for a^J, 
G/R m in is equal to 0.04 instead of 0.08 which would be in better 
agreement with our smallest values. 

6.4. Deep Impact 

The spectra of comet 9P/Tempel 1 were analysed before and af- 
ter the July 4 collision with the Deep Impact spacecraft. No dif- 
ference was observed in the intensity ratios and the lines widths 
before and after the impact (Fig. 15). We have no spectrum dur- 
ing the impact. The first spectrum after the impact was taken at 
53559.954 UT (i.e. ~6 hours after the impact). Cochran (2008) 
did not notice any change during the impact either. 
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Fig. 15. The G/R intensity ratio of comet 9P/Tempel 1 from MJD 
53553.955 to 53562.956 UT. The Deep Impact event corresponds to the 
vertical line (53559 UT). The average value is shown by the solid red 
line and the dashed lines correspond to the standard deviation (<x). The 
mean value around the impact (from 53553.955 UT to 53562.956 UT) 
is 0.04 ± 0.01. No variation of the ratio was noticed after the impact. 
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6.5. Solar activity 

Our data were obtained over the 23 th solar cycle characterized 
by its maximum of activity in 2001 and its minimum activity in 
2008. Despite the decrease of the solar activity during the obser- 
vations of C/2001 VI (NEAT) and 8P/Tuttle, no change in the 
intensities and widths of lines were noticed. Solar activity does 
not appear to have any influence on the results of our work. 



the low values of G/R obtained for comets 9P/Tempel 1 and 
8P/Tuttle. 
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7. Conclusion 

From January 2003 to September 2011, 12 comets of various 
origins were observed at different heliocentric distances and 48 
high resolution spectra were obtained with the UVES spectro- 
graph at VLT (ESO). Within this whole sample, we observed the 
three [OI] forbidden oxygen lines with high resolution and sig- 
nal to noise and we measured the intensity ratios (76300/^6364 an d 
G/R = /557 7 /(/6300 + ^364)) as well as the FWHM (v) i of the lines. 
The results can be summarized as follows : 

1. The /63oo/^6364 ratio (3.12 + 0.02) is in very good agreement 
with the branching ratio obtained by quantum mechanics, es- 
pecially the work done by Galavis et al. (1997). This result 
shows the high quality of our measurements. 

2. From theoretical values given in Table 1, the G/R ratio 
for comets observed below 2 AU (0.09 ± 0.02) confirms 
H2O as the main parent molecule photo-dissociating to pro- 
duce oxygen atoms. However, when the comet is located 
at larger heliocentric distances (r > 2.5 AU), the ratio in- 
creases rapidly showing that an other parent molecule is con- 
tributing. We have shown that CO2 is the best candidate. 
Measuring the G/R ratios could then be a new way to es- 
timate the abundances of CO2, a very difficult task from the 
ground. Assuming that only the photo-dissociation of H2O 
and CO2 produce [OI], we found a relation between CO2 
abundance and the heliocentric distance of the comet. The 
Q4 (NEAT) abundance of C0 2 at 3.7 AU is found to be 
-75%. 

3. The intrinsic green line width is wider than the red ones by 
about 1 km s _1 . Theoretical estimations considering Ly-or 
as the only excitation source for the two states lead to the 
conclusion that the excess energy for O('D) is larger than 
for O('S) which is in contradiction with our observations. 
This discrepancy might be explained by a different nature of 
the excitation source and/or a contribution of CO2 as parent 
molecule to the O('S) state. Indeed, Bhardwaj & Raghuram 
(2012) have shown that for the photo-dissociation of CO2, 
the main excitation source might rather be the 995-1165 A 
band. To check this hypothesis quantitatively, it would be 
necessary to estimate theoretically the excess energy for the 
oxygen atoms when the wavelength band is 995-1 165 A, ac- 
counting for the photo-dissociation of both H2O and CO2. 
The widths of the three [OI] lines are similar in Q4 (NEAT) 
at ~ 3.7 AU. This could be in agreement with CO2 being the 
main contributor for the three [OI] lines at large heliocentric 
distance. Other comets at large r have to be observed to test 
this hypothesis. 

4. More CO2 (and CO) abundance determinations, together 
with G/R oxygen ratios and line widths at different heliocen- 
tric distances, are clearly needed in order to give a general 
conclusion about the oxygen parent molecule. 

5. The C0 2 rich comet 103P/Hartley 2 (~ 20%) does not 
present a high G/R ratio normally expected from Eq. 12. We 
suggest a new value of 0.5% for that could also explain 
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